Design, optimization, and inference of multiphasic decay of infectious virus particles.

The loss of virus particles is typically considered to arise from a first-order kinetic process. Signals of deviations from this exponential decay are often de-prioritized. Here, we propose methods to evaluate if a design is adequate to evaluate evidence for multiphasic virus particle decay and to optimize the sampling times of decay experiments, accounting for uncertainties in viral kinetics. First, we evaluate 1500 synthetic scenarios of biphasic decays, with varying decay rates and initial proportions of subpopulations. Robust inference of multiphasic decay is more likely when the faster decaying subpopulation predominates insofar as early samples are taken to resolve the faster decay rate. Overall, we find that design optimization leads to a better precision of estimation while reducing the number of samples. It helps to estimate adequately the fastest decay in 54% of situations vs. 41% using a non-optimized design. We then apply these methods to infer multiple decay rates associated with the decay of ΦD9, an evolved isolate derived from phage Φ21. A pilot experiment confirmed that ΦD9 decay is multiphasic, but was unable to resolve the rate or proportion of the fast decay subpopulation(s). We then applied optimal design methods to propose new ΦD9 sampling times. Using this strategy, we were able to robustly estimate both decay rates and their respective subpopulations. Notably, we conclude that the vast majority (94%) of the population decays at a rate 16-fold higher than a slow decaying population. Altogether, these results provide methods to quantitatively estimate heterogeneity in viral decay.

Rapid bacteria-phage coevolution drives the emergence of multiscale networks.

Interactions between species catalyze the evolution of multiscale ecological networks, including both nested and modular elements that regulate the function of diverse communities. One common assumption is that such complex pattern formation requires spatial isolation or long evolutionary timescales. We show that multiscale network structure can evolve rapidly under simple ecological conditions without spatial structure. In just 21 days of laboratory coevolution, Escherichia coli and bacteriophage Φ21 coevolve and diversify to form elaborate cross-infection networks. By measuring ~10,000 phage-bacteria infections and testing the genetic basis of interactions, we identify the mechanisms that create each component of the multiscale pattern. Our results demonstrate how multiscale networks evolve in parasite-host systems, illustrating Darwin's idea that simple adaptive processes can generate entangled banks of ecological interactions.

Fecal virome transplantation is sufficient to alter fecal microbiota and drive lean and obese body phenotypes in mice.

The gastrointestinal microbiome plays a significant role in modulating numerous host processes, including metabolism. Prior studies show that when mice receive fecal transplants from obese donors on high-fat diets (HFD) (even when recipient mice are fed normal diets after transplantation), they develop obese phenotypes, demonstrating the prominent role that gut microbiota play in determining lean and obese phenotypes. While much of the credit has been given to gut bacteria, the impact of gut viruses on these phenotypes is understudied. To address this shortcoming, we gavaged mice with viromes isolated from donors fed HFD or normal chow over a 4-week study. By characterizing the gut bacterial biota via 16S rRNA amplicon sequencing and measuring mouse weights over time, we demonstrate that transplanted viruses affect the gut bacterial community, as well as weight gain/loss. Notably, mice fed chow but gavaged with HFD-derived viromes gained more weight than their counterparts receiving chow-derived viromes. The converse was also true: mice fed HFD but gavaged with chow-derived viromes gained less weight than their counterparts receiving HFD-derived viromes. Results were replicated in two independent experiments and phenotypic changes were accompanied by significant and identifiable differences in the fecal bacterial biota. Due to methodological limitations, we were unable to identify the specific bacterial strains responsible for respective phenotypic changes. This study confirms that virome-mediated perturbations can alter the fecal microbiome in vivo and indicates that such perturbations are sufficient to drive lean and obese phenotypes in mice.

Fecal virome transplantation is sufficient to alter fecal microbiota and drive lean and obese body phenotypes in mice.

Background: The gastrointestinal microbiome plays a significant role in numerous host processes and has an especially large impact on modulating the host metabolism. Prior studies have shown that when mice receive fecal transplants from obese donors that were fed high-fat diets (HFD) (even when recipient mice are fed normal diets after transplantation), they develop obese phenotypes. These studies demonstrate the prominent role that the gut microbiota play in determining lean and obese phenotypes. While much of the credit has been given to gut bacteria, studies have not measured the impact of gut viruses on these phenotypes. To address this shortcoming, we gavaged mice with viromes isolated from donors fed HFD or normal chow. By characterizing the mice’s gut bacterial biota and weight-gain phenotypes over time, we demonstrate that viruses can shape the gut bacterial community and affect weight gain or loss. Results: We gavaged mice longitudinally over 4 weeks while measuring their body weights and collecting fecal samples for 16S rRNA amplicon sequencing. We evaluated mice that were fed normal chow or high-fat diets, and gavaged each group with either chow-derived fecal viromes, HFD-derived fecal viromes, or phosphate buffered saline controls. We found a significant effect of gavage type, where mice fed chow but gavaged with HFD-derived viromes gained significantly more weight than their counterparts receiving chow-derived viromes. The converse was also true: mice fed HFD but gavaged with chow-derived viromes gained significantly less weight than their counterparts receiving HFD-derived viromes. These results were replicated in two separate experiments and the phenotypic changes were accompanied by significant and identifiable differences in the fecal bacterial biota. Notably, there were differences in Lachnospirales and Clostridia in mice fed chow but gavaged with HFD-derived fecal viromes, and in Peptostreptococcales, Oscillospirales, and Lachnospirales in mice fed HFD but gavaged with chow-derived fecal viromes. Due to methodological limitations, we were unable to identify specific bacterial species or strains that were responsible for respective phenotypic changes. Conclusions: This study confirms that virome-mediated perturbations can alter the fecal microbiome in an in vivo model and indicates that such perturbations are sufficient to drive lean and obese phenotypes in mice.

Comparison of bacterial suppression by phage cocktails, dual-receptor generalists, and coevolutionarily trained phages.

The evolution and spread of antibiotic-resistant bacteria have renewed interest in phage therapy, the use of bacterial viruses (phages) to combat bacterial infections. The delivery of phages in cocktails where constituent phages target different modalities (e.g., receptors) may improve treatment outcomes by making it more difficult for bacteria to evolve resistance. However, the multipartite nature of cocktails may lead to unintended evolutionary and ecological outcomes. Here, we compare a 2-phage cocktail with a largely unconsidered group of phages: generalists that can infect through multiple, independent receptors. We find that λ phage generalists and cocktails that target the same receptors (LamB and OmpF) suppress Escherichia coli similarly for ~2 days. Yet, a "trained" generalist phage, which previously adapted to its host via 28 days of coevolution, demonstrated superior suppression. To understand why the trained generalist was more effective, we measured the resistance of bacteria against each of our phages. We find that, when bacteria were assailed by two phages in the cocktail, they evolved mutations in manXYZ, a host inner-membrane transporter that λ uses to move its DNA across the periplasmic space and into the cell for infection. This provided cross-resistance against the cocktail and untrained generalist. However, these mutations were ineffective at blocking the trained generalist because, through coevolutionary training, it evolved to bypass manXYZ resistance. The trained generalist's past experiences in training make it exceedingly difficult for bacteria to evolve resistance, further demonstrating the utility of coevolutionary phage training for improving the therapeutic properties of phages.

Leapfrog dynamics in phage-bacteria coevolution revealed by joint analysis of cross-infection phenotypes and whole genome sequencing.

Viruses and their hosts can undergo coevolutionary arms races where hosts evolve increased resistance and viruses evolve counter-resistance. Given these arms race dynamics (ARD), both players are predicted to evolve along a single trajectory as more recently evolved genotypes replace their predecessors. By coupling phenotypic and genomic analyses of coevolving populations of bacteriophage λ and Escherichia coli, we find conflicting evidence for ARD. Virus-host infection phenotypes fit the ARD model, yet genomic analyses revealed fluctuating selection dynamics. Rather than coevolution unfolding along a single trajectory, cryptic genetic variation emerges and is maintained at low frequency for generations until it eventually supplants dominant lineages. These observations suggest a hybrid 'leapfrog' dynamic, revealing weaknesses in the predictive power of standard coevolutionary models. The findings shed light on the mechanisms that structure coevolving ecological networks and reveal the limits of using phenotypic or genomic data alone to differentiate coevolutionary dynamics.

Coevolutionary phage training leads to greater bacterial suppression and delays the evolution of phage resistance.

The evolution of antibiotic-resistant bacteria threatens to become the leading cause of worldwide mortality. This crisis has renewed interest in the practice of phage therapy. Yet, bacteria's capacity to evolve resistance may debilitate this therapy as well. To combat the evolution of phage resistance and improve treatment outcomes, many suggest leveraging phages' ability to counter resistance by evolving phages on target hosts before using them in therapy (phage training). We found that in vitro, λtrn, a phage trained for 28 d, suppressed bacteria ∼1,000-fold for three to eight times longer than its untrained ancestor. Prolonged suppression was due to a delay in the evolution of resistance caused by several factors. Mutations that confer resistance to λtrn are ∼100× less common, and while the target bacterium can evolve complete resistance to the untrained phage in a single step, multiple mutations are required to evolve complete resistance to λtrn. Mutations that confer resistance to λtrn are more costly than mutations for untrained phage resistance. Furthermore, when resistance does evolve, λtrn is better able to suppress these forms of resistance. One way that λtrn improved was through recombination with a gene in a defunct prophage in the host genome, which doubled phage fitness. This transfer of information from the host genome is an unexpected but highly efficient mode of training phage. Lastly, we found that many other independently trained λ phages were able to suppress bacterial populations, supporting the important role training could play during phage therapeutic development.

Fecal Viral Community Responses to High-Fat Diet in Mice.

Alterations in diet can have significant impact on the host, with high-fat diet (HFD) leading to obesity, diabetes, and inflammation of the gut. Although membership and abundances in gut bacterial communities are strongly influenced by diet, substantially less is known about how viral communities respond to dietary changes. Examining fecal contents of mice as the mice were transitioned from normal chow to HFD, we found significant changes in the relative abundances and the diversity in the gut of bacteria and their viruses. Alpha diversity of the bacterial community was significantly diminished in response to the diet change but did not change significantly in the viral community. However, the diet shift significantly impacted the beta diversity in both the bacterial and viral communities. There was a significant shift away from the relatively abundant Siphoviridae accompanied by increases in bacteriophages from the Microviridae family. The proportion of identified bacteriophage structural genes significantly decreased after the transition to HFD, with a conserved loss of integrase genes in all four experimental groups. In total, this study provides evidence for substantial changes in the intestinal virome disproportionate to bacterial changes, and with alterations in putative viral lifestyles related to chromosomal integration as a result of shift to HFD.IMPORTANCE Prior studies have shown that high-fat diet (HFD) can have profound effects on the gastrointestinal (GI) tract microbiome and also demonstrate that bacteria in the GI tract can affect metabolism and lean/obese phenotypes. We investigated whether the composition of viral communities that also inhabit the GI tract are affected by shifts from normal to HFD. We found significant and reproducible shifts in the content of GI tract viromes after the transition to HFD. The differences observed in virome community membership and their associated gene content suggest that these altered viral communities are populated by viruses that are more virulent toward their host bacteria. Because HFD also are associated with significant shifts in GI tract bacterial communities, we believe that the shifts in the viral community may serve to drive the changes that occur in associated bacterial communities.

De novo origins of multicellularity in response to predation.

The transition from unicellular to multicellular life was one of a few major events in the history of life that created new opportunities for more complex biological systems to evolve. Predation is hypothesized as one selective pressure that may have driven the evolution of multicellularity. Here we show that de novo origins of simple multicellularity can evolve in response to predation. We subjected outcrossed populations of the unicellular green alga Chlamydomonas reinhardtii to selection by the filter-feeding predator Paramecium tetraurelia. Two of five experimental populations evolved multicellular structures not observed in unselected control populations within ~750 asexual generations. Considerable variation exists in the evolved multicellular life cycles, with both cell number and propagule size varying among isolates. Survival assays show that evolved multicellular traits provide effective protection against predation. These results support the hypothesis that selection imposed by predators may have played a role in some origins of multicellularity.